"Centrifugal fragmentation" in the photodissociation of in intense laser fields 
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By means of quantum-dynamical and classical trajectory calculations of photodissociation in 
strong laser fields, it is shown that for certain combinations of pulse durations and intensities the 
rotational dynamics can lead to centrifugal fragmentation. In that case, the photofragments exhibit 
characteristic angular distributions. The classical calculations provide a transparent physical picture 
of this mechanism which is also very well established in collisions between atomic nuclei or liquid 
droplets: non-rotating systems are stable, whereas rotating systems fragment due to the decrease 
of the fragmentation barrier with increasing angular momentum. 



Centrifugal fragmentation (CF) is a universal phe- 
nomenon of rotating masses. It happens if the angu- 
lar momentum of the rotating object exceeds a critical 
value l CI and the system breaks into two pieces due to 
centrifugal forces. In the seminal work of Cohen, Plasil 
and Swiatecki [l[ , these critical i-values have been calcu- 
lated for nuclei, droplets and gravitating masses within 
the universal rotating liquid drop model (RLDM). Later, 
critical angular momenta for metallic clusters have been 
estimated, too For nuclei and droplets, the criti- 
cal angular momenta can be created and measured in 
heavy ion cofiisionsj3[ and in collisions between macro- 
scopic droplets (j4a|, respectively. During these colli- 
sions, an intermediate, fast rotating compound is formed 
which undergoes fission for I > l CI and remains stable for 
I < l CI . The ^ cr -values are obtained from the universal 
-^-dependence of the fusion cross sections as function of 
the incident energy E (see dill for details). The exper- 
imental Z cr -values for nuclei and droplets are in excellent 
agreement with the predictions of the RLDM, although 
their absolute values differ by more than 20 orders of 
mag nitude (0 and references therein). 

Diatomic molecules exposed to strong laser pulses may 
also rotate because the interaction of the induced or per- 
manent dipole moment of the molecule with the laser 
field leads to an effective torque towards the laser polar- 
ization axis. To what extent rotational effects or even 
CF may influence the photodissociation process, how- 
ever, is a challenging, up to date problem (71-lllj and the 
central topic of this work. This will be illustrated for 
the photodissociation dynamics of nature's most simple 
molecule, Hj, where alignment @j[l2-14| plays an im- 
portant role for the interpretation of the experimentally 
observed angular distributions of the fragments (see, e. g., 
the summary given in ref. 



11]) 



Molecular alignment is often conveniently discussed in 
terms of light-dressed, adiabatic Floquet-surfaces (l5l - fl9| 
as a function of the internuclear distance R and the angle 
9 between the molecular axis and the laser polarization 
axis. One commonly distinguishes between two different 



alignment mechanisms: 

(i) Geometric alignment. Upon excitation from the at- 
tractive Is cr g ground state to the repulsive Is a u state, 
a fragmentation channel opens at 9 = 0° leading to pref- 
erential dissociation of molecular ions initially aligned 
along the polarization axis via the bond softening (BS) 
mechanism on the lower Floquet surface. In contrast, the 
upper surface exhibits a barrier at 9 = 0° (bond hard- 
ening (BH) mechanism (USUI) and a minimum at 
9 = 90° leading to "counter-intuitive" angular distribu- 
tions [H. 

(ii) Dynamic alignment. The laser field generates an ef- 
fective torque towards the laser polarization axis which 
leads to rotation of initially non-aligned molecules and 
results in photofragments near 9 « 0° via BS on the 
lower Floquet surface [H [H [H-IH . 

However, experimentally, it is difficult to distinguish be- 
tween both mechanisms, in particular, because they oc- 
cur simultaneously with fragments at 9 ~ 0°. 

In this work, it will be shown that, under certain 
conditions, the rotational dynamics can lead to an ad- 
ditional photodissociation mechanism. By solving the 
nuclear time-dependent Schrddinger equation (TDSE) a 
pronounced and characteristic double humped angular 
distribution of the fragments is found. The first peak 
centered along the polarization axis results from the well 
known BS mechanism of geometrically and dynamically 
aligned molecules. The second peak, however, deviating 
15° — 20° from the polarization axis, can decisively be 
attributed to a "centrifugal fragmentation" mechanism 
where the molecule accumulates sufficient angular mo- 
mentum that the centrifugal force drives the fragments 
over the dissociation barrier. Under suitable conditions 
the CF-peak even dominates over the BS-contribution at 
9 = 0°. 

The dynamics of in an intense laser field e(t), in 
dipole approximation and length gauge, is described by 
the time-dependent Schrddinger equation (TDSE) for the 
electron coordinate (r) and the internuclear distance vec- 
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tor (R) 

i-|*(R, r, t) = [f R + H e (r; R) + &(<)] tf (R, r, t) (1) 

where Tr is the kinetic energy of the nuclei, -ff e (r;R) 
the electronic Hamiltonian. The electric field e(t) = 
eof(t) cosiut is linearly polarized along the z-axis and 
characterized by the amplitude eo> constant carrier fre- 
quency w and the envelope function f(t). For this we use 
a sin 2 -shaped turn-on of 15 fs to simulate a cw-laser and 
f(t) = sin 2 ^ for a pulsed laser with total duration T. 
Atomic units (a.u.) are used unless specified otherwise. 

We expand the full wavefunction ^(R, r, t) in paramet- 
rically R-dependent eigenstates ^{(r; R) of the electronic 
Hamiltonian as 



*(R,r,t) = ^i^(R,t)$ £ (r;R 



(2) 



with the pulse length T. For the cut-off Rd, we have cho- 
sen values of 4 a.u. and 10 a.u. for low and high vi- 
brational states, respectively. They lie beyond the BS 
dissociation barrier and have been chosen such that fur- 
ther increase does not affect the distribution Pjj. 

The angular momentum distribution of the fragments 
as a function of time Pi(i), i.e., the weights of the par- 
tial waves I contained in the wavepacket, is obtained by 
projecting the nuclear wave function fij on spherical har- 
monics Y™ 

oo 

R d «' m 
7T 2n 2 

x / dQ sine* J #F, ro *(6l,^)0 ? (R,t;m) . (6) 



taking into account the two lowest, i.e., the a g and a u 
states in the actual calculations, and thus neglect ioniza- 
tion. We have checked that the inclusion of the higher ix u 
state does not change the results for the laser parameters 
used in this paper. 

Within the framework of the Born-Oppenheimer (BO) 
expansion in eq. ([2]), the nuclear time-dependent 
Schrodingcr equation reads 



ij^(R,i) 



i d 2 



L 



2[idR 2 2nR 2 



vdR) 



n£(R,t) 

(3) 



with the eigenvalues of the electronic Hamiltonian 
V^(R) and the dipole matrix elements D^ n (R) = 
J d 3 r$^(r; R)z$ l; (r; R). The initial state is chosen as 
a product of the vibrational eigenstate v, and spherical 
harmonics with fixed angular momentum U , starting from 
the CTg-surface 

n ag {R,6^;p,l i ,m;t = 0)= X u,i i (R)-YlT( e ^) • ( 4 ) 

Equation ([3]) is then solved numerically for all m-values 
— li < m < k using the code WavePacket [23| by 
expanding the wave function in a basis of plane waves 
and spherical harmonics, and applying the split opera- 
tor scheme 28(. For the radial grid, an equally-spaced 
grid with 4096 points ranging from 0.2 a.u. to 236 a.u. 
is used. Angular momenta up to I = 59 are taken into 
account. The time step is 10 -2 fs. 

The angular distributions of the photofragments Pd (0) 
are calculated from the nuclear wave function according 
to 



Pd(9) = 



1 



2L 



dR 



fi ? (R,T;m) 



(5) 



In addition to the BO-expansion @ we solve the TDSE 
(HJ) by expanding the total wave function *P(R, r, t) in 
time-dependent Floquet states [Io| and consciously take 
into account only the lowest Floquet surface Vp(R, t) 
(adiabatic Floquet approximation [29|). The results are 
compared to those of the full BO-calculations. In doing 
so we are able to distinguish BS and CF from other disso- 
ciation mechanisms (BH, multiphoton dissociation). In 
this approximation, the nuclear wave function Qp(R, t) 
is determined by a single TDSE 
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n F (R,t) . 

(7) 
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FIG. 1: Fragment angular distributions Pd(0) (left, solid 
lines), and corresponding angular momentum distributions 
Pi(t) (right, linear density plots) of dissociating H^, ini- 
tially in its rotational and vibrational ground state, in a 
cw-laser with A = 230 nm, I = 70 TW/cm 2 (a,b) and 
I = 180 TW/cm 2 (c,d). The dotted lines in (a), (c) cor- 
respond to the reduced one-Floquet-surface calculations (see 
text). 
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The angle- and time-dependence of the Floquet surface 
Vf(R, t) is given implicitly by an effective field strength 
e eS (9,t) = e f{t)cos9 as V F (R,t) = V F {R,e oS (9,t)). 
The numerical integration of (|7|) is performed in a similar 
way as for ©. 

First, we consider the wavepacket dynamics of H^" ini- 
tially in its ground state {v = 0, = 0), exposed to a 
cw laser with a wave length of A = 230 nm and two in- 
tensities / = 70 TW/cm 2 and / = 180 TW/cm 2 (Fig. 
[1]). To ensure the cw-character, i.e., convergence in cal- 
culated probabilities, a propagation time of T = 300 
fs is sufficient. At / = 70 TW/cm 2 a forward peaked 
angular distribution is obtained as expected from the 
ordinary BS-mechanism on the lowest Floquet-surfacc 
with geometric and dynamic alignment contributions. At 
/ = 180 TW/cm 2 , however, a second peak at angles be- 
tween 9 w 15°... 20° occurs. The time-dependence of the 
angular momentum distributions in Fig. Q] (b) and (d) 
suggests that the appearance of this peak is directly re- 
lated to distinctly higher I- values. To exclude any other 
mechanism, caused by multiphoton dissociation or BH ef- 
fects, we compare the results of the full BO-calculations 
with those obtained in the ground-state Floquet approxi- 
mation ([7]) . Evidently, the results of both calculations are 
identical and, thus, centrifugal effects are unambiguously 
the origin of the second maximum in Pd{9) ("centrifugal 
fragmentation" ) . 

Second, we employ a pulsed laser keeping wavelength 
and intensity unchanged (Fig. [2]). For short and long 
pulses (T = 30 fs and T — 240 fs) single-peaked angu- 
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lar distributions are obtained. In the fast growing elec- 
tric field of the short pulse (Fig. [2] (a)) dissociation of 
non-aligned molecules becomes also possible resulting in 
a broad distribution. On the other hand, the slow in- 
crease of the field strength in long pulses (Fig. [2] (e)) 
supports strongly the dynamical alignment effect which 
obviously results in a very narrow distribution with a 
distinct larger maximum value at 9 = 0°. For intermedi- 
ate pulse lengths (T = 85 fs), however, the bimodal dis- 
tribution is recovered, associated with distinctly larger 
angular momenta of the fragments as compared to the 
other cases (Figs. [2] (b),(d),(f)). The perfect agreement 
between the exact and restricted one-Floquet surface cal- 
culations confirms again the one photon nature of the 
dynamics in all three cases. Obviously, for the interme- 
diate pulse length the molecule is rotationally accelerated 
to overcome the dissociation barrier on the first Floquet 
surface at 9 — 0°; the fragments, however, are scattered 
to final angles 9 ss 15. ..20° due to their large angular 
momenta (see also the discussion in terms of classical 
trajectories, Fig. |4]). 

Whereas for in its rotational and vibrational 
ground state (v = 0, k = 0) the CF effect shows up as a 
relatively small contribution in the total angular distribu- 
tion Pd{9) (note the logarithmic scale in Figs. [T] and [2]), 
it can dominate the fragmentation mechanism for state- 
selected molecules (h > 0) where li = 2 is the leading 
component in the thermal distribution for para-H^. This 
is demonstrated in the upper part of Fig. [3] where, for 
fixed laser parameters, the angular distributions for 
in its rotational and vibrational ground state (Fig. [3] (a)) 
and excited (Fig. E](b)) Hj" are compared. The strongly 
forward peaked distribution for U = characterizes the 




FIG. 2: Same as in Fig. [TJ but for a pulsed laser with 
T = 30 fs (a,b), 85 fs (c,d), 240 fs (e,f) and A = 230 nm, 
/ = 180 TW/cm 2 . 



FIG. 3: Typical shapes of the angular distributions Pd{9) 
belonging to different dissociation mechanisms (BS, BH and 
CF) emerging in the fragmentation of depending on the 
initial state v, U (solid lines). 

The dotted lines correspond to the reduced one-Floquet- 
surface calculations. The Frank-Condon-averaged distribu- 
tion is shown as well (lower right). The laser parameters are 
A = 266 nm, T = 130 fs and I = 120 TW/cm 2 . 
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ordinary BS mechanism (Fig. [3] (a)). Somewhat surpris- 
ingly, a minimal shift in the initial angular momentum 
(li = 2) leads to a qualitative change in the shape of the 
distribution with a distinct maximum at finite angles as 
the fingerprint of the CF mechanism (Fig. [3] (b)). It is 
important to note that the CF maximum in Fig. [3J (b) 
disappears for longer (300 fs) and shorter (30 fs) pulse 
durations. This behavior demonstrates again the dynam- 
ical nature of the effect (see also Fig. [2]). 

In addition, we consider vibrationally excited 
molecules (Fig. [3J (c),(d)). For the state selected case 
with v = 4, li = 2 a very broad wavelike distribution with 
an abrupt decrease for angles 9 > 60° is observed (Fig. [3] 

(c) ). The complete disagreement with the isotropic dis- 
tribution, predicted by the ground-state Floquet calcu- 
lation (dotted line), proves definitely that this angular 
distribution is due to the dominating influence of the 
excited Floquet surface, i.e., the "counterintuitive" BH 
fragmentation mechanism [20|. Hence, the occurrence, 
and thus the experimental observability of the CF mech- 
anism, depends also on the initial ro-vibrational state of 
the molecule. 

In many experiments, the neutral H2 molecule serves 
as the precursor for the fragmenting Hj" and, thus, a dis- 
tribution of different z^- values weighted by their Franck- 
Condon factors [3(| is examined at once. For this case, 
the predicted angular distribution is shown in (Fig. [3] 

(d) ). As can be seen for the actual laser parameters, 
the CF mechanism "survives" the averaging procedure 
and leads to an angular distribution with a dominating 
CF-maximum at finite angles. At the same time, the BH 
mechanism generates the long tail contribution for angles 
0>3O° (cf. Fig. 1(b), (d)). 

Complementary to the full (|3]) and restricted 
wavepacket dynamics we have also performed purely clas- 
sical calculations where the nuclei are propagated on 
the lowest time-dependent Floquet surface Vp with the 
Hamiltonian 



p2 r2 

2^1 2fiR 2 



V F (R,t) 



(8) 



The resulting Hamiltonian equations of motion for the 
momenta (Pr,L) and coordinates (R,9) represent the 
classical analogue to Schrodinger equation ([7]). 

Two representative trajectories (i?,6>) are shown in fig- 
ure |4] An initially almost aligned molecule (with initial 
conditions Qi = 5°, Ri = 1.8 a.u.) essentially vibrates 
and slightly librates within the potential and remains sta- 
ble. A non-aligned molecule with 9i — 60°, Ri — 2.4 a.u., 
however, vibrates and strongly rotates simultaneously, 
accumulating a large angular momentum and overcom- 
ing the barrier at the saddle point near 9 = 0°, R « 3 a.u. 
to finally fragment toward an angle of 9 w —15°. Also af- 
ter averaging over trajectory bundles, the CF mechanism 
is recovered, but overestimated with regard to the quan- 
tum calculations due to the lack of tunneling. However, 




1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 
R(a.u.) 

FIG. 4: Classical picture of centrifugal fragmentation. Up- 
per: Floquet potential surface Vf{R, 9) of with two sam- 
ple trajectories for A = 230 nm, / = 70 TW/cm 2 and 
T — 300 fs (cw). Lower: The effective potential V B s{R) = 



V F (R)- 



at 6 = 0° for angular momenta I — and I — 9. 



The fragmentation barrier is clearly smaller for I — 9 due to 
centrifugal forces. For comparison, the energy of the lowest 
vibrational state is plotted as straight line. 



these calculations provide a clear and simplified phys- 
ical picture of the underlying mechanism of CF: Non- 
rotating (aligned) molecules are stable, whereas rotating 
(hence initially non-aligned) molecules fragment owing to 
the decrease of the fragmentation barrier with increas- 
ing angular momentum (compare the effective potentials 
V e e(R) for I = and I = 9 in the lower part of Fig. HJ). 

Exactly the same mechanism leads to the instability 
(i.e. CF) of fast rotating nuclei, clusters and droplets [l|, 
2j, where the critical angular momenta are described by 
the RLDM. In the present molecular case, the ^ cr -values 
depend on the laser parameters via Vp (R, 9) and their 
systematic investigation remains an interesting topic for 
future theoretical studies. 

In summary, by means of quantum and classical cal- 
culations, we have definitely shown, that the rotational 
dynamics of H^ in intense laser fields can lead to centrifu- 
gal fragmentation. The effect shows up in characteristic 
side-peaked angular distributions of the fragments, pro- 
vided optimal laser parameters are chosen to accelerate 
the molecule rotationally above a critical Z-value. Exper- 
imentally, the observability of CF may be complicated by 
uncertainties in intensity and focal volume of the laser as 
well as finite temperature effects, and thus represents a 
great challenge. 

This research is funded by the Danish National Re- 
search Foundation's Center for Molecular Movies and by 
DFG through program SFB 450. 
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